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Abstract 

We calculate the production rate of neutral fermions in linear magnetic fields through the Pauli 
interaction. It is found that the production rate is exponentially decreasing function with respect 
to the inverse of the magnetic field gradient, which shows the non-perturbative characteristics 
analogous to the Schwinger process. It turns out that the production rate density depends on both 
the gradient and the strength of magnetic fields in 3+1 dimension. It is quite different from the 
result in 2+1 dimension, where the production rate depends only on the gradient of the magnetic 
fields, not on the strength of the magnetic fields. It is also found that the production of neutral 
fermions through the Pauli interaction is a magnetic effect whereas the production of charged 
particles through minimal coupling is an electric effect. 
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I. INTRODUCTION 



It is well known that the interaction of charged spin-1/2 fermions with the electromag- 
netic field is described by the minimal coupling in the form of Dirac equation. Pauli [l| 
suggested a non-minimal coupling of spin 1/2-particle with electromagnetic fields, which 
can be interpreted as an effective interaction of fermion to describe the anomalous mag- 
netic moment of fermions. Pauli interaction is particularly interesting for describing the 
interaction of neutral particle with electromagnetic field provided it has a non-vanishing 
magnetic dipole moment ^. One of the immediate possibilities |4j is the electromagnetic 
interaction(Pauli interaction) of neutrinos, which are recently confirmed to have a non-zero 
mass with mixing jjj. The presence of the magnetic dipole moment implies that neutrino 
can directly couple to the electromagnetic field, which leads to a variety of new processes 

DO 

One of the interesting phenomena with the strong electromagnetic field configuration is 
the pair creation of particles. The well known example is the Schwinger process with the 
minimal coupling, in which charged particles are created in pairs 7] under a strong electric 
field. However it has been demonstrated that no particle creation is possible under the 
pure magnetic field configuration even with the spatial inhomogeneity j^. For a neutral 
particle with the Pauli interaction, the inhomogeneity of the magnetic field coupled directly 
to the magnetic dipole moment plays an interesting role analogous to the electric field for a 
charged particle. The non-zero gradient of the magnetic field can exert a force on a magnetic 
dipole moment such that the neutral fermion can get an energy out of the magnetic field. 
Hence it will affects the vacuum structure greatly for the strong enough magnetic field 
as for the case of charged particles in the strong electric field. It is then interesting to 
see whether the vacuum production of neutral fermion with a non-zero magnetic moment 
in an inhomogeneous magnetic field is possible on the analogy of the Schwinger process. 
Interestingly it has been demonstrated in 2+1 dimension that the magnetic dipole coupled 
to the field gradient induces pair creations in a vacuum jjj. In this work, we will present a 
realistic calculation in 3+1 dimension for the pair production rate of neutral fermions with 
non- vanishing magnetic dipole moment. 
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II. PRODUCTION RATE OF NEUTRAL FERMIONS THROUGH PAULI IN- 
TERACTION 



The simplest Lagrangian for the neutral fermion, which couples to the external electro- 
magnetic field, was suggested by Pauli long time ago The Dirac equation with Pauli 
term is given by 

C = $V+£o*'F ia ,-m)i/>, (1) 

where a^ u = | py M , j u ], g^ u = (+, — , — , — ). fi in the Pauli term measures the magnitude of 
the magnetic dipole moment of fermion. Pauli term can be considered as an effective inter- 
action which describes anomalous magnetic moment of fermion or as an effective magnetic 
moment induced by the bulk fermions in a theory with large extra dimensions Q|. The 
corresponding Pauli Hamiltonian operator is 

H = a ■ (p — ifj,/3E) + f3(m — fia ■ B), (2) 

where a 1 = }^e % ^ k a^ k . If the magnetic field is stronger than the critical field, B c = m/fi, a 
negative energy state of m — \fiB\ appears. Thus, we will consider magnetic fields weaker 
than the critical magnetic field, B < B c . 

In general, the effective potential, V e s(A), for a background electromagnetic vector po- 
tential, Afj_, can be obtained by integrating out the fermion: 

% I d 4 xV eS (A[x}) = Trln{^ + ^F tiV -m)-^—}, (3) 

where F^ = d^A u — d u A^. The decay probability of the background magnetic field into the 
neutral fermions is related to the imaginary part of the effective potential V e g(A), 

p — 1 _ | e */ d 4 xV cfi (A[x])^ _ ^ _ e -2/m / d 3 xdtV cS (A[x]) 

That is, the twice of the imaginary part of the effective potential V^Ajx]) is the fermion 



production rate per unit volume 



ll|: w(x) = 2Im{V eS {A[x})). 



Using the charge conjugation matrix C: 

C^C'^-rf, Ca^C- l = -a T »\ (5) 



and the identity 



b Jo s 



we can write the effective potential Vrff(A[a;]) as follows 

V eS (A[x}) = - — e - ism2 tr{< x\e is ^ + ^ F ^ 2 \x > - < x\e isp2 \x >). (7) 
2 Jo s 

For an inhomogeneous magnetic field, we consider a static magnetic field configuration 
of z-direction with a constant gradient along ^-direction, B = B(x)z, such that 

F u = B(x) = B + B'x = B'x, (x = x + x, x = —). (8) 

B 

It is not necessary to consider an infinitely extended ever- increasing magnetic field to meet 
the linear magnetic field configuration. Because the particle production rate density is a 
local quantity, it is sufficient to have a uniform gradient magnetic field in the Compton 
wavelength scale of the particle. 

For this background magnetic field, we get 

(J + ^F, u ) 2 = -pi -pj + (jjlB'x + ifp + *7°P 3 ) 2 - , (9) 

where 7 3 = 7°7 1 7 2 and 7 = 7 1 7 2 7 3 . 

The second term in Eq.((7j) is easily calculated to be 

=tr < x\e lsp2 \x >= (10) 

Inserting a complete set of momentum eigenstates, the first term of Eq.(JJ|) can be written 
by 



/ A) = tr < x \e~ is{p2+p2 ~^ B '* +i ^ P0+i ^° P3)2+ ^ B ' y2} 
1 ,7T , 1 



X > 



is J 

e TSr (7°P3+7 3 P0 ) (Pi -Pi ) e -is^B'j 2 ^ ^ 

where we have used [7 s , 7 2 ] = = [7 , 7 2 ], [pi, x] = —i and 

Pi + v\ ~ {n-B'x + ry 3 po + «7°P3) 2 = 

e -^(7°P3+7 3 P0)pi el x pi( p 2 +p 2 _ ( AtjB ' x )2) e ^^°P3+7 3 P0)pi e - J ,. 0Pl> (12) 

Using the properties of gamma matrices 

{ 7 ~0 ; 7 ~3} = 0, tr(7 3 7 2 ) = = tr(7° 7 2 ), (13) 



we get 



tre ^^°P3+7 3 P o)K-Pi) e - t ^^ = 4 CQsh(g/ij g/) CQS { (Pi Pi) (pg _ p 2 )1/2} (M) 

/iir 

Thus, f ^ can be written as follows 



(A) 



4 , 7T , 1 



(2vr) 4 



cosh(s//5') / dp 1 dp' 1 dp dp 3 e i{pi - p '^ x 



cos {^B^(Po -^) 1/2 } < Pi|e-^ W Vi > . (15) 

The matrix elements in momentum space in Eq. (jl5)l correspond to the matrix elements 
of the evolution operator for the simple harmonic oscillator with an imaginary frequency, 
oj = 1i\iB' and m—\. Then we get 

< pxle-^-^'^V! >= J dx'dx" < Pl \x" > U(x", s; x', 0) < x'\p\ >, (16) 

where U(x, s; x', 0) is given by 



U(x, s; x', 0) ee< x, S | e - Js(p2+ ^ * V , >= ( ) 1/2 e luJ - . (17) 

4m sm out 

Performing the x' and x" integration explicitly, we get 

< Pi e iPl M Pi >= ( — ) 2 e 8a ( 1 - cos ^ s )e 8c "( 1+c ° s " s ) , (18) 

' 7r 2asinc<js 

where a = 



4 sin ujs ' 

Inserting Eq. (jl8|) into Eq.()15|l. the integration over pi and p[ gives 



v 



(A) 



2 . 7T . 1 



(2tt) 4V 2s 

where, defining 7, a and p_ as follow 



)5 C oth(s/i J B')/C, (19) 



coth(su-B') 1 . 9 9 . 1 , . . 

7^ 4 ^ % o=— (pg-pg)s, p^ft-pi, (20) 



/C is given by 

K = ~ r dp.dp dp 3 e iS:p -(e iap - +e- iap ~)e 
2 J— 00 

1 7T 1 /o /* ■ (a+a) 2 . (g-a) 2 

-f—W 2 / dp dp 3 (e ^ +e ^ 



2 v ry' 

47r(^') 2 (^)^ - 2( M 5) 2 m(-)^ /' d£(l - £)e* 2 . (21) 
z «7 jo 
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T? {i jf * (1 ~^' pBfe " - s + ^k*™ 2 -- (24) 



Thus, -i/ A ) is reduced to the following 

= — L[±{(s/iB') coth^E')} 1 + 2(/i£) 2 -{(s y uE') coth(s/i5')}5 /* (1 - 0e^ ?2 ] 

47T 2 S 2 S JO 

(22) 

and finally we get the effective potential V e g as follows 

V eS = Z - r^e—\v^-v^). (23) 

j£ j s 

The effective potential for the uniform field configuration can be obtained by putting 
\iB' = to get 

The divergent contributions at s = are removed by adding local counter terms of (fiB) 2 , 
and (fiB) 4 . This implies the renormalization of the magnetic moment \i to the measured 
value and the coupling of (/i-B) 4 to zero presumably The effective potential, Eq.pijl. for 
uniform magnetic fields is found to be real, which implies a stable magnetic background 
of w = 0. For magnetic fields weaker than the critical field B c = m/fi, using a contour 
integration in the fourth quadrant, the integration can be done along the negative imaginary 
axis giving the finite real effective action as 

n, = r % { I + _ /■ , m - (25) 

4n z Jo s z 2 12 Jo 
The leading radiative correction term for a weak B field is 



5V eS = {fi f 2 . (26) 
240n 2 m 2 



For an inhomogeneous field configuration, fiB' ^ 0, the effective potential is given by 
V* = f ${V^-oth(^ S ) / de(l - 0e^-*"> - 1 + ^}e"^ 



+ 3^/ o °° J{(^coth(^)) 3 / 2 - 1 - l^|£l! } e— (27) 



8tt 2 

where an additional divergent contribution at s = is removed by adding a local counter 
term of (/i-B') 2 m the second term. 

The leading radiative correction terms of the effective potential for a small gradient weak 
B field are calculated as given by 

cS 240vr 2 m 2 2887r 2 m 4 487r 2 m 2 960vr 2 m 4 ' 1 J 
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It is found that the effective potential, Eq.(|27|). has a non- vanishing imaginary part, 
which implies that the background of inhomogeneous magnetic field configuration is un- 
stable against the creation of neutral fermions with Pauli interaction. From the imagi- 
nary part of the effective potential Eq. (J27|) . we obtain the production rate density, w(x) = 
2Im(V eS (A[x})). 

2 2 

Introducing dimensionless parameters defined as v = sfiB', A = ■n^pr, K = , the 
production rate density w(x) in the unit of the fermion mass is finally given by 

2m 4 f°° dv A 1 Xv 

w(x) = — — — — / — {v v cothvF(— tanht>, Xv) cos Aw sinAv} 

Att 2 XkJo v 2X k k ' ; 2 12k ; 

r^{(,; CO tht;) 3 / 2 -l--}sinAt;, (29) 
Jo v 6 2 



4tt 2 A 2 



where 



F(a,b) = f 1 dm - cosK 2 - b) 
Jo 

= {sin(a — b) + sin(6)} 

2a 

+ J— {cos(6) FresnelC(W— ) + sin(6) FresnelS(W — )}. (30) 
V 2a \j n \ n 

Since the scale of inhomogeneity less than Compton wavelength of the fermion is irrelevant 
to the particle production through this process, we take in this work the spatial gradient of 
the magnetic field \B'\ to be smaller than the ratio of field strength, \B\, to the Compton 
wavelength — , that is, A > y/n. 

The integration Eq. (j29j) is a finite integration, but has singularities along the imaginary 
v axis similarly to the Schwinger's result, where the residue calculation gives the analytic 
WKB type expression. However, the integration Eq. ()29j) has essential sigularities along the 
imaginary axis, so that it seems not possible to get a usual analytic WKB type expression 
using a contour integration. Therefore we use numerical integrations to investigate the 
properties of the production rate density given by Eq.(|29p. For the numerical calculation, 
we consider the case of k > 1 and A > 2 as an example. Numerical integrations of the 
production rate density show that the second term of Eq.(|29jl is negligible compared to the 
first term for A > 2. The production rate shows exponential monotonic decrease for A > I^Jk 
and k > 1. The production rate, w(x), is calculated for k — 1.0 and 2.0 as a function of A. 
The results in the unit of are as shown in [Fig. 1-2]. The results of numerical calculations 
are represented by dots in the figures. The numerical integrations of Eq. (j29|) suffer from 



large oscillatory fluctuation, which is an unavoidable feature due to the violent oscillations 
in the integrand as discussed in jl^]. To get an analytic expression, we obtain the best fit of 
the numerical results to the curves in the form of s^e~ bKX for [Fig. 1-2]. We can observe that 
the particle creation rate is an exponentially decreasing function with respect to the inverse 
of the field gradient, 

w ^ e -constantxm 2 /|/i.B'|^ 

which shows the characteristics of the non-perturbative process. This can be understood 
as a quantum tunnelling through a potential barrier of height ~ 2m of a particle exposed 
to an potential energy ~ fi\B'\x due to the inhomogeneous magnetic field coupled to the 
magnetic dipole moment through Pauli interaction. It is similar to the Schwinger process 
of electron-positron pair creation in the strong electric field, where the creation rate is 
decreasing exponentially [ill, w ~ e ~ coristantxm2 /\ eE \_ 

One can see that the production rate is suppressed very rapidly when the field strength 
becomes weaker than ~ m/fi as well as the inhomogeneity scale is bigger than Compton 
wave length scale. For the inhomogeneity in the Compton wave length scale, the rate per 

4 

unit time per unit volume is typically of order ^ for the critical field strength. 

However for the realistic estimation of the production rate, more precise information 
on the mass and magnetic moment of a particle and the strength of the magnetic field as 
well as the scale of the spatial inhomogeneity of the field in consideration are needed for 
the observational possibility. As a possible environment, let us consider the pair creation 
of neutrino with a non-zero magnetic dipole moment in the vicinity of the very strongly 
magnetized compact objects with B = 10 15 G as a typical strength Q|. Taking the 
possible magnetic moment to be as large as the experimental upper bound fi v = 10 _11 //b and 
the mass of the neutrino to be m u ~ 10~ 2 eV constrained by the solar neutrino observations, 
the critical field strength is estimated to be B c ~ 10 17 G. One can see that the condition for 
k > 1 or B < B c assumed in this work is satisfied. Since the the scale of the inhomogeneity 
is naturally about the size of the compact object, R ~ 10 4 m, which is much larger than 
the Compton wave length ~ lCT 4 m, the production rate is expected to be substantially 
suppressed from the typical rate such that it may not provide sufficiently high luminosity 
for the neutrino detectors. 
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III. DISCUSSION 



We have examined the vacuum production of neutral fermions in inhomogeneous mag- 
netic fields through Pauli interaction. The fermions, which are coupled to the background 
electromagnetic field through Pauli interaction, are integrated out and there appears an 
imaginary part in the effective action. It turns out that the production rate density depends 
on both the gradient and the strength of magnetic fields in 3+1 dimension, which is quite 
different from the result in 2+1 dimension [9], where the production rate depends only on 
the gradient of the magnetic fields, not on the strength of the magnetic fields. The difference 
can be attributed to the different nature of spinors in 3+1 and 2+1 dimensions. 

The vacuum production of fermion with the Pauli interaction is found to be a magnetic 
effect. Explicit calculations with a linear electric field configuration of x— direction with 
a constant gradient along x— direction such that E x = E + E'x shows that the effective 
potential has no imaginary part when the singularities are regularized properly. It can be also 
shown by substituting B —>■ iE and B' — > iE' in the effective potential for a pure magnetic 
field Eq. (l23j) with the s integration along the imaginary axis. Therefore one can see that 
the pair creation through Pauli interaction is a purely magnetic effect. It is an interesting 
result when compared to the pair creation of charged particles through the minimal coupling, 
which is known to be an electric effect jlli ] . 

Although the production rate density in this work has been derived for fiB' = constant, 
it can be applicable to various types of magnetic fields provided that the magnetic field is 
linear in the scale of Compton wavelength of the particle considered because the particle 
production rate density is a local quantity. It may be therefore applicable to a spatially 
slowly varying fiB'(x) as a good approximation if the gradient variation is very small in 
the Compton wavelength scale. For the realistic calculation of the production rate, more 
precise information on the mass and magnetic moment of a particle and the strength of the 
magnetic field as well as the scale of the spatial inhomogeneity of the field in consideration 
are needed for the observational possibility. 
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